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ABSTRACT 

A novel way of looking at the evolution of star clusters is presented. With a dynamical 
temperature, given by the mean kinetic energy of the cluster stars, and a dynamical 
luminosity, which is defined as the kinetic energy of the stars leaving the cluster in 
analogy to the energy of photons emitted by a star, the dissolution of star clusters is 
studied using a new dynamical temperature-luminosity diagram for star clusters. The 
investigation contains a parameter-space study of open clusters of up to TV = 32768 
single-mass stars with different initial density distributions, half-mass radii, tidal con- 
ditions and binary fractions. The clusters show a strong correlation between dynamical 
temperature and dynamical luminosity and most of the investigated cluster families 
share a common sequence in such a dynamical temperature-luminosity diagram. Devi- 
ations from this sequence are analyzed and discussed. After core collapse, the position 
of a cluster within this diagram can be defined by three parameters: the mass, the 
tidal conditions and the binary fraction. Due to core collapse all initial conditions are 
lost and the remaining stars adjust to the given tidal conditions. Binaries as internal 
energy sources influence this adjustment. A further finding concerns the Lagrange radii 
of star clusters: Throughout the investigated parameter space nearly all clusters show 
a constant half-mass radius for the time after core collapse until dissolution. Further- 
more, the ratio of half-mass radius to tidal radius evolves onto a common sequence 
which only depends on the mass left in the cluster. 

Key words: stellar dynamics - galaxies: star clusters - methods: ./V-body simulations 
- Hertzsprung-Russel (HR) diagram - open clusters and assoziations: general 
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1 INTRODUCTION 

In 1913 Henry Norris Russell presented his work on a rela- 
tion between the spectral classes of stars and their absolute 
magnitude at a meeting of the Royal Astronomical Society 
l|Russelllll913l ). The diagram he showed later became famous 
as the Hertzsprung-Russel Diagram (Ejnar Hertzsprung was 
the first to anticipate the existence of a relation between the 
two quantities) and became one of the most important tools 
for the study and understanding of stellar evolution. 

The temperature- luminosity diagram, which is a deriva- 
tive of the original Hertzsprung-Russel diagram, shows a 
tight relation between a star's temperature and its luminos- 
ity for the main-sequence phase of stars. 

A star cluster is a system which shows some analogies 
to a star. First of all, the stars within a cluster follow a 
velocity distribution, which is established through two-body 
encounters, just like a gas or plasma does through collisions 
of particles. Therefore it is possible to assign a dynamical 
temperature, T, to the stars in a cluster (Sec.|3j. Secondly, 



a star cluster constantly loses a certain fraction of its stars 
through escape, like a star constantly emits photons. It thus 
appears plausible to define a dynamical luminosity, L, in 
terms of the energy carried away by the stars per unit time 
(Sec. . 

To find if a dynamical temperature-luminosity relation 
exists for dynamical systems like star clusters is the motiva- 
tion of this work, because it is likely that it would prove very 
useful for describing global cluster properties and evolution. 
Using numerical simulations, the behaviour of star clusters 
within a dynamical T-L diagram (Sec. [5]) and the influence 
of different initial conditions on its development (Sec. [6} is 
studied. 

Throughout the performed parameter-space study the 
half-mass and tidal radius of the clusters are investigated 
in detail, as they turned out to not behave as expected. 
The half-mass radius stays constant for a large fraction of a 
cluster's life-time and the ratio of half-mass radius to tidal 
radius evolves along a single sequence, independent of initial 
conditions. 
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But first of all a reference model, the "Standard Clus- 
ter" , is defined and investigated in detail (Sec. [2}, which will 
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help establishing a dynamical temperature and a dynamical 
luminosity. 



2 STANDARD CLUSTER 

Since the number of possible initial models for N-body com- 
putations of star clusters is infinite, it is necessary to scan 
the available parameter space systematically. Therefore it is 
highly convenient to define a "Standard Cluster" with care- 
fully chosen, plausible initial conditions, which is taken as a 
reference model. After detailed studies of this type of model, 
the effects of changes to the initial conditions can be traced 
efficiently. From these studies of the different variations, gen- 
eral features of cluster evolution can be deduced. 

For this investigation about one hundred and fifty clus- 
ters with different initial conditions were calculat ed. They 
all were com puted with the iV-body code NBODY4 |Aarsethl 
1 19991 . 120031 ) until total dissolution (i.e. until N = 10). Al- 
though NBODY4 is capable of including stellar evolution 
and similar effects, they were not implemented here to focus 
on pure dynamical evolution. 

It turned out that for clusters with up to ~2000 stars, 
a regular workstation is sufficiently fast. The larger models 
were computed on the GRAP E-6A special purpose comput- 
ers of the Argelander Institute (|Fukushige. Makino fc Kawail 
120051 '). 

The Standard Cluster is defined as follows: 

(i) Plummer radial density profile, 

(ii) half-mass radius of 0.8 pc, 

(iii) 1000 equal-mass stars of 1 Mq , 

(iv) tidal field, 

(v) no primordial binaries. 

Details to the several choices are given below. 

(i) The Pl ummer profile is chosen because of its analytical 
conve nience l|Aarseth. Henon fe Wielen 1974; He ggie fc Hutl 
l2003h . 

(ii) T he choice of R h — 0.8 pc is motivated by iKroupal 
(1995a|), who finds by inverse dynamical population synthe- 
sis that the majority of Galactic field stars is probably born 
in aggregates of this size. 

(iii) N = 1000 is a good compromise between statistics 
and computing time; on a regular workstation such a cluster 
can be computed within a few hours. 

Constraining the masses of the stars to one solar mass puts 
the focus on the pure dynamical evolution of the clusters, 
since stellar evolution needs not to be treated. Furthermore 
it avoids further parameters which have to be fixed initially 
(e.g. the initial mass function). 

(iv) The tidal field is added in the near field approx- 
imation, which is a linearisat ion of the galaxy potential 
l|Ross. Mennim fc He ggie 1997j). Here it is assumed that the 
cluster moves on a circular orbit in the Galaxy, which is 
treated as a point mass M G = 9.565 • 10 10 M at dis- 
tance Rg = 8.5 kpc, corresponding to a rotational speed 
of Vrot = 220kms~ 1 . The Standard Cluster is therefore sub- 
ject to boundary conditions typical for clusters in the solar 
neighbourhood. 

The computations are conducted in a rotating reference 



frame, which moves with the cluster's initial angular veloc- 
ity, with the origin of the coordinate system lying in the 
initial barycentre of the cluster. The rc-axis points to the 
Galactic centre and the i/-axis is parallel to v ro t, while the 
z-axis is p erpendicular to the orbit of t he cluster. For further 
details see iFukushige fc Heggie! (|2000h . 

All models start in dynamical equilibrium but the tidal 
field is added after setting up the initial positions and veloc- 
ities of the stars. Thus a few stars may already be unbound 
at the beginning. Because the initial clusters are quite small 
compared to the tidal radius this effect may be neglected, 
but becomes more important for initially very extended clus- 
ters. 

(v) Binaries slow down computations enormously because 
time step s have to be ver y small for their integration. Hence, 
although IKroupal (|l995ar ) suggests a rather large primordial 
binary fraction, the models did not contain primordial bina- 
ries, i.e. the primordial binary fraction, fbin, given by 

fbin = N^+W s ' (1) 
where Nun is the number of binary systems and N s is the 
number of single stars, such that N = 2Nbin + Ns, was set 
to zero. 

Since the Standard Cluster offers little statistics due to 
its small N, accurate results on statistical properties can 
only be found by combining many renditions of the same 
model, differing only in the initial positions and velocities 
of the stars. Therefore in the following sections the insights 
from simulations of 27 renditions of the Standard Cluster 
are presented, eventually supported by exemplary outcomes 
from other models. In Sec. |6] each variation is discussed in 
detail. 



3 DYNAMICAL TEMPERATURE 

A powerful analogon to the dynamics of a cluster with 
stars of a single mass is giv en by the kinetic theory of 
gases (|Louis fc Spurzemlll99ll ). In this framework the clus- 
ter is treated as an ideal gas with point-particles interacting 
only through collisions. Although this analogy exhibits some 
flaws, since interactions between stars are not confined to the 
short moment of collision, one can get valuable quantitative 
results out of it. 

In analogy to an ideal gas, a dynamical temperature of 
a cluster can be defined through the mean velocity of the 
member stars. The three-dimensional squared velocity dis- 
persion is given by 

o- 2 = — ^ ((v xi — v x ) 2 + (vyi — v y ) 2 + (v zi — Vz) 2 ) , (2) 

i 

where N is the number of bound stars (which will be defined 
in Sec. [4}, the v's are the velocity components and the v's 
the corresponding mean velocities. For a star cluster at rest, 

— 1 — 1 2 

Ekin = -mv 2 = -ma , (3) 

where m is the mass of the stars, which is taken to be the 
same for all stars, and Ekin is their mean kinetic energy. 
Furthermore, it is possible to define the pressure of a star 
cluster, 
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Figure 1. Upper panel: Lagrange radii of a Standard Cluster; 
LR2, LR10 and LR90 are respectively the 2%, 10% and 90% La- 
grange radii (the vertical line shows t cc , the horizontal line is the 
fitted mean value of the half-mass radius, i?^, for t > t cc ). While 
the core undergoes oscillations, the outer shells expand until they 
reach a nearly constant state (especially the half-mass radius). 
The tidal radius (dotted curve) decreases continuously due to on- 
going mass loss. The middle panel shows the mass evolution, and 
the bottom panel the temperature evolution with and without 
replacing the binaries with their centre-of-mass particles. After 
rising slightly, the temperature drops after core collapse. Without 
replacement, binary systems dominate the cluster temperature. 
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Figure 2. Temperature evolution of a cluster with the same prop- 
erties like a Standard Cluster, except fn n = 0.95, plotted with 
and without correcting for binaries. Tight binaries completely 
dominate the cluster temperature when they are not replaced 
by the corresponding centre-of-mass particles. 
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2V 



(4) 



where JV is the number of particles and V the volume of 
the system. With the so found quantity P and the ideal gas 
equation, 



PV = NkT, 



(•>) 



where k is Boltzmann's constant, a temperature T can be 
defined, which is proportional to the squared velocity dis- 
persion of the cluster, i.e. the mean kinetic energy of the 
member stars: 



1 2 

-ma 



kT. 



(6) 



In the case of a star cluster the temperature is usually de- 
noted as kT, thus has the dimension of an energy. Like for 
a perfect gas, the temperature therefore gives a direct idea 
of the distribution of velocities within the cluster. 

But this approximation breaks down when two parti- 
cles form a binary, since then the premise of an ideal gas 
is not valid any more. A single tight binary is able to com- 
pletely dominate the cluster temperature, due to the large 
internal stellar motion, and distort the frame of kinetic the- 
ory. Therefore these systems have to be neglected for the 
measurement of the temperature. 

Taking a look at the actual binary content of a Standard 
Cluster reveals a rather small fraction of cluster stars cap- 
tured in binaries, i.e. the number of binaries in a Standard 
Cluster is mostly less than ten which validates the above 
approximations (see also Fig. [4] in Sec. [4j . Hence, for the 
temperature analysis, i.e. in eq. [5J binaries are replaced by 
their centre-of-mass particles. 

This corrected temperature is a direct measure of the 
velocity distribution of the stars within the cluster. This dis- 
tribution is supposed to be Maxwellian, which implies that 
there is a high-velocity tail of stars, and, since a star cluster 
has a finite escape speed, that a number of stars are un- 
bound. When these hot stars have left the cluster, the left- 
over stars will reestablish a new velocity distribution with a 
lower temperature within a two-body relaxation time. This 
process continues until the cluster has dissolved. This can 
be seen in the lower panel of Fig. [T] where the value of kT 
of a Standard Cluster in the course of time is shown with 
and without replacing the binaries. In this context the im- 
portance of the temperature correction becomes clear when 
looking at the temperature evolution of a cluster with a 
substantial fraction of primordial binaries (Fig. [5]). With- 
out replacing the binaries by their centre-of-mass particles, 
kinetic theory completely breaks down, since the measured 
temperature is not related to the actual velocity distribution 
of the stars within the cluster. In this particular example, 
the temperature of the cluster at dissolution time would be 
even higher than in the beginning, which is to be expected 
due to the negative heat capacity of gravity. 

The temperature evolution of the Standard Cluster can 
be explained by looking at the top and middle panel of 
Fig. [T] In the beginning the cluster is very compact and 
sits deeply within the tidal sphere. As the core immediately 
starts collapsing (top panel) the temperature rises and the 
outer layers expand. Core collapse is reached at the indi- 
cated time, t cc , where the Standard Cluster has expanded 
sufficiently to fill its tidal radius and starts spilling over the 
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Figure 4. Upper panel: velocities of escaping stars as a function 
of time. While the majority of escapers leave the cluster with a 
velocity comparable to the escape velocity and Eki n /kT < 5.1 
("Evaporaters"), a few stars are virtually shot out of the cluster 
("Ejecters"). The dotted vertical line gives the time of core col- 
lapse, t cc , which coincides with the occurrence of the first ejecter. 
Lower panel: smoothed number of bound binary systems. Dur- 
ing core collapse, binaries are formed which then start producing 
ejecters. 



tidal radius. This can also be observed by the change of slope 
of the mass evolution (middle panel). 

A further feature in the top panel of Fig. [T] is striking: 
after core collapse the half-mass radius is nearly constant 
for the rest of the cluster's life-time. This feature can be ob- 
served in almost all computed clusters and will be discussed 
in detail in Sec. [6] The clusters adopt to the tidal condi- 
tions and settle to an equilibrium half-mass radius, where 
energy production in the core (which causes cluster expan- 
sion) and energy loss at the tidal radius balance each other. 
The dashed horizontal line in Fig. [T] gives a fit to the half- 
mass radius after core collapse until dissolution time, tdis , of 
this particular cluster. This has been done for all computed 
models. As can be seen in Tab. Q] (column R[), all clusters 
with N = 1000 and Roai = 8.5 kpc establish a half-mass 
radius of about 2 pc. 



4 DYNAMICAL LUMINOSITY 

In general the luminosity of a system can be defined as the 
amount of energy emitted in a given time interval. For a star 
this energy is given by the number of emitted photons and 
their energy. 

When trying to define a dynamical luminosity for star 
clusters in terms of kinetic energy carried away by escaping 
stars, it is first necessary to define when a star has left the 
cluster. Answering this question is not trivial since stars can 
still be bound even beyond the tidal radius or when they 



alrea dy have the necessary energy for escape (I Baumgardt 
2001). A sufficient and simple escape criterion is (jAarseth 

booa i 



x d \ > 2R tide , 



(7) 



where Xi is the star's position, Xd is the cluster's density cen- 
tre and Rude is t he tidal radius. The latter can be calculated 
as (|Spitzerll 19871 ) 



Ri 



Rg 



f M \ 1/3 
\3M G ) 



(8) 



Hence, the easiest way of measuring a cluster's luminos- 
ity is to put a virtual sphere around the cluster with radius 2 
Rude an d record the kinetic energy of each star leaving this 
sphere within a given time interval. In this way the luminos- 
ity is determined on the "surface" of the cluster. This also 
implies that all quantities, like the temperature, are most 
appropriately determined for all stars within this sphere. 

Investigating the velocities of t he escaping stars reveals 
two distinct families of escaper s (G lcrsz & Heggid Il994bl; 
iBaumgardt, Hut fc Heggi3l2002ft . In the left panel of Fig. [3] 
the distribution of the kinetic energies of the escaping stars 
of all computed Standard Clusters is shown. Here the ener- 
gies have been normalised by the current temperature kT of 
the cluster, which is reasonable to do when trying to com- 
pare different kinds and states of clusters, because in a hotter 
system the mean velocity is shifted to higher values and so 
is the mean escape velocity. In this plot n esc is normalised 
such that 



= 1. 



(9) 



Two lognormal distributions which overlap at about 
Ekin/kT — 5.1 are readily apparent. The origin of these 
two families of escapers lie in the way they are produced. 

(i) The slow escape is usually called evaporation. It is due 
to a number of weak encounters, which gradually increase 
the star's energy until it reaches a velocity higher than the 
escape velocity. Evaporation therefore is a consequence of 
the high-velocity tail of the velocity distribution of stars 
within a cluster, which is continuously established through 
two-body relaxation. Since, for a given cluster, the relax- 
ation time decreases with increasing temperature and vice 
versa, it is supposed that the evaporation rate is correlated 
to the temperature. 

(ii) High-velocity escapers are due to ejection processes. 
These can be encounters between a star and a binary or two 
binaries, where a binary system is hardened while the excess 
energy is transferred to the other star(s). This implies that 
the energy taken away by ejecters is mostly stored in a re- 
maining binary system, unless this binary recoils sufficiently 
strongly to leave the cluster. As the stars have to come very 
close to each other, these events are quite rare and mostly 
happen in the cluster core. 

For further considerations it is convenient to split esca- 
pers into ejecters and evaporaters. For the Standard Clus- 
ter a threshold {Eki n /kT) limit can be determined from the 
intersection point of the two lognormal distributions. Of 
course, this also has to be checked for clusters with differ- 
ent intital conditions. For all computed classes of models 
normalised velocity distributions were acquired, all resem- 
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Figure 3. Distribution of mean normalised escape velocities of 27 Standard Clusters (left) and 9 isolated clusters (right). n esc gives the 
probability of detecting an escaper with the given Ekin/kT. The dotted lines show the fitted lognormal distributions which intersect at 
about Ekin/kT = 5.1. In isolated clusters evaporation is less common than in Standard Clusters. 



bling the distribution of the Standard Cluster: Here also 
two lognormal distributions were fitted and their intersec- 
tion point was determined. The individual results can be 
found in Tab. □ 

Since all models show a limit of about the same value 
like the Standard Cluster, a weighted mean of these values is 
suggested as the threshold between evaporaters and ejecters, 

(TF) =5.1 ±0.2. (10) 

V n>£ / limit 

The same analysis has even been done for the computed iso- 
lated clusters, where no tide re moves the slow stars, h ence 
evaporation is less efficient (see iBaumgardt et al.1120021 ). In 
Fig. [3] the relative increase of ejecters, in number and also 
in relative energy, can be seen. The fit yields a limit of 
Ekin/kT = 4.3 ± 0.6 which agrees within the errorbars with 
the results for clusters in tidal fields and therefore seems to 
be independent of the tidal conditions. 

Looking at the velocities of escaping stars reveals that 
most ejections take place at and right after core collapse 
(see Fig. |4}, which is highly correlated to the formation of 
binaries. The probability of a close encounter is highest in 
this phase. 

In fact, the first ejecter of a star cluster can be taken as 
a trigger for core collapse, i.e. just at the moment when the 
first escaper with Ekin/kT > 5.1 is detected, the core has 
reached its densest phase. This criterion should be preferred, 
because it is less arbitrary than the usual way of defining 
t cc , which is the formation of the first hard binary with a 
given amount of binding energy dGiersz fc Heggiei [l994bl ). 
The first ejecter on the other hand marks a definite change 
in the core energy evolution, since a large amount of energy 
is removed immediately. 

The considerations made so far suggest the definition 
of a dynamical total luminosity, L to t, which is given by the 
kinetic energy carried away by escaping stars. Since two dis- 
tinct processes are at work it is reasonable to split this total 



luminosity into an evaporation luminosity, L ev , and an ejec- 
tion luminosity, L e j , such that 

Ltot — L ev + L e j. (11) 

The evaporation luminosity of a cluster is comparable to 
the luminosity of a star, where photons carry away heat 
with energies proportional to their frequency. So, to stay in 
the framework of stars, evaporaters may be the analogon 
to black-body radiation, while ejecters axe more like solar 
flares. Evaporation takes place constantly and the amount of 
energy being lost by it strongly depends on the cluster's tem- 
perature, much like the luminosity of a main-sequence star 
correlates with its effective temperature. Ejection is mostly 
due to structural effects happening deep inside the cluster, 
such as core collapse, thus due to the formation and hard- 
ening of binaries. It therefore happens more randomly, but 
is a good tracer for binary activity. 

But a last paradoxon still has to be solved: if an escap- 
ing star carries away positive energy and the total energy 
is negative in the beginning, the cluster should get increas- 
ingly bound in the course of time. How come that clusters 
dissolve finally? The answer is that the binding energy is 
redistributed among the left-over stars predominantly such 
that binaries are formed or hardened. These binaries store 
binding energy until they are removed from the cluster due 
to a three- or more-body event in the cluster core. In Fig. [5] 
the energy evolution of a Standard Cluster can be seen. 
The jumps indicate binary-escaper events, which all coincide 
with a high- velocity ejecter of more than 10 km/s going in 
the opposite direction of the binary-escaper. This supports 
the argument that ejecters are strongly correlated to binary 
activity. As can be seen in the figure, by simply counting the 
steps in the energy evolution, this particular Standard Clus- 
ter produces 12 binary escapers within its life-time (com- 
pared to about 40 ejecters and more than 900 evaporaters) . 
Hence binary escapers are too rare to define a separate lu- 
minosity. So, if a binary escapes, just the kinetic energy of 
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Figure 5. Upper panel: Total energy of the cluster with (thin 
dashed curve) and without (thick solid curve) counting the inter- 
nal energy of binaries (binaries being replaced by centre-of-mass 
particles). Lower panel: total energy per star. From being tightly 
bound in the beginning, the cluster dissolves completely. While 
stars constantly leave the cluster the binding energy is predom- 
inantly stored in binary systems. The jumps in energy coincide 
with binary escapers which take away large fractions of the total 
binding energy. The dashed vertical line marks core collapse. 



its centre of mass is measured and added to the appropriate 
luminosity. 



5 THE TEMPERATURE-LUMINOSITY 
DIAGRAM 

Having now introduced the dynamical temperature, kT, and 
the dynamical luminosity, L, a Hertzsprung-Russel-type di- 
agram for star clusters can be drawn (Fig. [6} . 

Taking only escapers with E^in/kT < 5.1 into account 
and averaging over all 27 Standard Clusters yields a smooth 
curve in the dynamical temperature-luminosity diagram. The 
Standard Cluster starts deep within its tidal radius with an 
evaporation luminosity of about 0.1 M0km 2 s _2 Myr _1 and a 
temperature of ~ 2M0km 2 s -2 , then increases its luminosity 
before core collapse, after which it cools down constantly 
until it has dissolved completely to a temperature of about 
zero. 

The total luminosity, including evaporation and ejec- 
tion, is a factor of 10 higher and is governed by single 
ejecters. It therefore does not yield a smooth curve but shows 
a variety of peaks which are of statistical nature. For more 
massive clusters, which show pronounced core oscillations, 
these peaks should be correlated to these oscillation, because 
during contraction phases the central density, and hence the 
probability of few-body encounters, is highest. 
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Figure 6. Mean temperature-luminosity diagram of all 27 Stan- 
dard Clusters taking only evaporation into account (L ev ) and 
taking both, evaporation and ejection, into account (Ltot)- Er- 
rorbars represent the standard deviation from the mean. Since 
ejection only takes place in dense environments, it is dominant in 
the beginning but decreases as the cluster dissolves. 



Since in iV-body models with single-mass particles the 
primary parameter for defining the state of the cluster is the 
total mass rather than the age, the different renditions are 
binned, combined and averaged at equal masses in steps of 
50 Mq. The errorbars give the standard deviation of single 
models from the mean. As is evident, the ejection luminosity 
suffers from much stronger statistical fluctuations than the 
evaporation luminosity. 

Finally, which deductions can be drawn from the com- 
putations of the Standard Cluster? First of all, a reasonable 
way of defining a dynamical temperature and a dynamical 
luminosity have been found. This allows to set up a dy- 
namical temperature-luminosity diagram and look at cluster 
evolution in a completely new way. The individual Standard 
Clusters evolve along a single sequence within this diagram, 
which is in this case more comparable to a cooling sequence 
of white dwarfs than to a main sequence. This is due to the 
cluster temperature, which is constantly declining, except 
for the period before core collapse. 

Furthermore, as can be seen in the different slopes of the 
two curves in Fig.[6j energy loss through ejection is most im- 
portant around core collapse and then decreases compared 
to evaporation. This difference is due to a different depen- 
dence of the two luminosities on cluster parameters such as 
the number of stars and the density. While the ejection lu- 
minosity is mainly related to the central velocity dispersion 
and density, because binaries are most active in the core, the 
evaporation luminosity depends on the total density distri- 
bution, because evaporaters are generated everywhere in the 
cluster IjBaumgardt et al.l 2002). 

After this pioneering work, the variations of the ini- 
tial parameters and their influences on the cluster evolution 
within the dynamical temperature-luminosity diagram can 
be traced. 
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6 VARIATIONS 

The following variations of the Standard Cluster were inves- 
tigated (the full list of model properties is given in Tab.[l|: 

(i) As an alternative to the Plummer profile, several mod- 
els with an initial King profile (Wo = 5) were computed 
l|Kindll966h . These models were chosen to have the same 
initial half-mass radius like the Standard Cluster of 0.8 pc, 
hence they are tidally underfilling by a factor of three. 

(ii) The half-mass radius was varied to 0.4, 1.6, 2.4, 3.2 
and 4.0 pc. 

(iii) The number of stars was varied from 500 up to 32768 
stars, corresponding to 500-32768 Mg. The effect of adding 
a mass function has to be studied in some future work, but 
is expected to have a large influence on cluster evolution (se e 
lHenonlll969l ; Ivan Albadalll968l ; Ide La Fuente Marcos! fl995h . 

(iv) Models were calculated at a larger Galactic radius of 
Rg = 85 kpc. For these models the tidal radius is about 60 
pc, which is five times larger than for a Standard Cluster. 
The cluster therefore can expand much further which, at the 
same time, increases its dynamical time scale significantly. 

Furthermore isolated models were computed to get a bet- 
ter idea of the general influence of a tidal field. 

(v) To see the effect of a substantial binary fraction, N- 
body integrations with ft, in = 0.95 were performed. No other 
binary fractions were investigated to keep the number of 
models containing primordial binaries to a minimum be- 
cause of the long computation times for such models, which 
are up to 100 times larger than for models without primor- 
dial binaries. 

To minimise the changes from model to model with regard 
to future experiments, where a realistic multi-mass spectrum 
will be implemented, the binding energy distribution, fip, 
was specified in such a way that the binding energies of 
the binary systems follow the observed distribution of real 
binaries, but here applied to single-mass systems. 

For a realistic binary system, take two stars mi and m,2 
independently out of an init ial canonical mass distribution 
£(m) oc m~ ai l|Kroupall200"ll ) with 



qi = 1.3 for 0.08 sC m/M < 0.5, 
q 2 = 2.3 for 0.5 < m/Mg < 150. 



(12) 
(13) 



Chose a binary per iod P which foll ows the period distribu- 
tion fip derived bv lKroupal (|l995ch . 



ftp = 2.5 



IP - 1 



45 + (ZP-l) 



with 1 sC IP «; 8.43, 



(14) 



where IP = log 10 (P) and P is given in days in this equation. 
The semi-major axis, a, becomes 



a G[m\ + m-i) 

Finally, the binding energy of the binary is 



Ebin — 



2a 



(15) 



(16) 



The distribution of binding energies is then given by a triple 
integral over the mass ranges of the two composite stars and 
over the above given range of periods, 

fiE = J J J fiP i(m{]£(mz)8[lP -IP) d mi dm 2 d/P, (17) 

7711 77T2 IP 

with the substitution 



• 'v!H log 



(m 1 m 2 ) a 
mi + mi 



-IE, 

2 ' 



(18) 



where IE = log w (\E),i n \). The result is a distribution func- 
tion for binding energies, fiE, combining an initial mass 
function £(m) and a period distribution function fip. It is 
normalised such that 



fw(lE) dlE - 



(19) 



For the equal-mass binary population used in the present 
work a binding energy is chosen from the energy distribu- 
tion function fiE (Fig. [7J. The eccentricities, e, are then 
generated from the thermal eccentricity distrib ution func- 
tion, /(e) = 2e (e.g. iDuquennov fc Mavoilll99ll ). 

6.1 Density profile 

Fig.|S]depicts that, as expected, the tidally underfilling King 
models do not differ significantly from the Standard Cluster. 
Within the first few crossing times the initial density profile 
is completely lost, the King models therefore follow the same 
cooling track as the Standard Cluster. Furthermore, the half- 
mass radius of those models settles to the same equilibrium 
value of about 2 pc (see Tab.[T]). 

6.2 Half-mass radius 

Due to virial equilibrium, reducing the initial half-mass ra- 
dius raises the initial temperature of the system. A model 
with half-mass radius of 0.4 pc therefore exhibits a temper- 
ature which is twice the temperature of a Standard Cluster, 
while a model with 1.6 pc starts with half of this tempera- 
ture. 

But since these clusters have the same tidal radius, they 
redistribute their mass through two-body relaxation such 
that after core collapse the initially smaller or larger clus- 
ters do not differ in their properties (e.g. temperature, den- 
sity profile, binary fraction) from a Standard Cluster of the 
corresponding mass, except for their age. This holds espe- 
cially for the half-mass radius (Fig. [9}, which settles to the 
same equilibrium value of 2 pc as the Standard Cluster (see 
Tab.[T}. This indicates that the evolution after core collapse 
is independent of the initial density profile, which can also 
be nicely seen in the dynamical T-L diagram: after starting 
to "shine" they all share a common cooling track (Fig. I10[) . 

A cluster with an initial half-mass radius larger than 
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Figure 8. Mean temperature-luminosity diagram of clusters with 
an initial King density profile (Wo = 5, tidally underfilling by 
a factor of three). Upper panel: Evaporation + ejection. Lower 
panel: Evaporation. For comparison the thick grey lines give the 
cooling track of the Standard Cluster (Fig. IS) . 



the equilibrium value of 2 pc, e.g. 3.2 pc, however, evolves 
very differently to a Standard Cluster (Fig. [TTJ . Since the 
initial median two-body relaxation time, given by 

N 1/2 R V2 



Uel = 0.138— -j- ± , (20) 

M 1/2 G 1 / 2 InA 

where m is the mean ste llar mass and InA is the Coulomb 
logarithm (Spitzer 1987 j3, is about ten times larger than 
for a Standard Cluster, evolution towards core collapse will 
take much longer. In fact, the ratio of core collapse time 
and intital relaxation time is supposed to be a constant 
with value about 17 for Plummer models l|Baumgardt et all 
2002), hence it takes 10 times longer in such models to reach 
core collapse, while the dissolution time is barely changed 
(compare with Tab. [p. During this time the half-mass radius 
will be shrinking until the equilibrium value is established. 




1 iGiersz fe~ Heggic (1994a) find that the Coulomb logarithm can 
be approximated by InA ~ ln(0.11A r ). 



1000 900 800 700 600 500 400 300 200 100 
M[MJ 

Figure 9. Half-mass radii of models with initial of 0.4, 0.8, 
1.6, 2.4, 3.2 and 4.0 pc. All converge to an equilibrium value of 2 
pc after core collapse, when the clusters have lost about 30-50% 
of their initial mass. 



This has important consequences for the evolution of the 
other cluster parameters. 

The middle panel of Fig. [TT] shows, that mass loss is 
more rapid at the very beginning, because the larger cluster 
already fills the tidal sphere. This effect becomes more im- 
portant the larger the cluster is. But, as energy production 
and hence expansion is not so efficient, the Standard Cluster 
dissolves faster in the end. 

As in all other models, no binary ejecters are produced 
before core collapse. The total energy of these models is 
therefore nearly constant over a large time span, so is the 
temperature (lower panel of Fig. Hip . 

This has significant influence on the evolution in the dy- 
namical temperature-luminosity diagram. Since the temper- 
ature stays constant for about half of the cluster's lifetime, 
it spends much time at the same position on the common 
cooling track, until its structure has reached the Standard 
Cluster configuration and starts moving along with it un- 
til dissolution (Fig. I12|) . This is an excellent proof for the 
existence of a tight relation between temperature and lu- 
minosity. Furthermore this shows that the cooling sequence 
is also a dynamical main sequence, since star clusters with 
half-mass radii larger than the equilibrium value stay at the 
same spot on the sequence until their internal configuration 
has adopted to the tidal conditions, which can take up to 
half of their total life-time. 

The clusters with an initial Rh of 4.0 pc lose so many of 
their stars in the very beginning due to their initial dimen- 
sion, that their luminosity is initially slightly above the dy- 
namical main sequence. More extended clusters would have 
a even higher initial luminosity before they finally reach the 
dynamical main sequence. 

6.3 Number of stars 

Increasing the initial mass of the cluster causes a deeper po- 
tential well and hence a higher temperature. However, all 
clusters evolve along a common cooling sequence indepen- 
dently of their initial masses (Fig. I13p . In Fig. Q3] the re- 
sults of the most massive clusters are shown. This was done 
seperately for clarity reasons because for N = 32768 and 
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Figure 10. Mean temperature-luminosity diagram of clusters 
with initial half-mass radii of 0.4 and 1.6 pc. Upper panel: Evapo- 
ration + ejection. Lower panel: Evaporation. For comparison the 
thick grey lines give the cooling track of the Standard Cluster 
(Fig. U . 



N = 16384 only one model each was computed, which re- 
sults in a larger scatter. Nevertheless the two clusters clearly 
extend the cooling sequence to higher temperatures. Investi- 
gating the luminosity in dependence of the mass shows that 
after core collapse all models evolve on a common track, 
hence are fully characterised by the mass left in the cluster 

(Fi g . Hsu. 

This is also emphasized in Fig. 1161 which shows the ratio 
of half-mass radius to tidal radius as a function of the mass 
left in the cluster. It can be seen that all models start tidally 
underfilling and increase their Rh/ Rude initially, after which 
they evolve asymptotically onto a common sequence, inde- 
pendent of their starting condition. On this sequence the 
ratio Rh/Rtide increases, which could be interpreted in the 
way that, as the half-mass relaxation time (eq. [20} decreases 
with decreasing mass, energy production within the half- 
mass radius gets more efficient, leading to a nearly constant 
equilibrium half-mass radius opposed to a declining tidal 
radius. 
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Figure 11. Above: Evolution of the Lagrangian radii for a clus- 
ter with initial = 3.2 pc. Core collapse takes about half the 
cluster's whole lifetime. The middle panel shows the mass evolu- 
tion of the same cluster compared to a Standard Cluster, which 
is dissolving more slowly at the very beginning. In the bottom 
panel the temperature of the cluster is shown. During evolution 
towards core collapse the temperature stays nearly constant. 



This supports the assumption that, after core collapse, 
the state of a cluster in a given tidal condition just depends 
on its mass. All initial conditions are lost through core col- 
lapse and a universal density distribution is established. 

This finding of a constant Rh shows that the simplified 
theory of cluster evo l ution set-up for illustrative purposes by 
iBinnev fc Tremamel (|l987l , p. 525) does not capture cluster 
evolution correctly. According to this ansatz the half-mass 
radius decreases with ongoing mass loss, as 



Rh(t) — Rh 



M(t) 
M 



(21) 



where and Mo are the initial values of the half-mass ra- 
dius and the cluster mass. In this ansatz it is assumed, that 
evaporation is the dominant mode of escape which leaves 
the energy of the cluster nearly constant, as evaporaters take 
away almost no energy. This holds true for initially very ex- 
tended clusters for the time before core collapse, since no 
ejecters have been produced yet. If the cluster is too con- 
centrated in the beginning, i.e. the initial half-mass radius 
is smaller than the equilibrium value, the cluster expands 
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Figure 12. Mean temperature-luminosity diagram of clusters 
with initial half-mass radii of 2.4, 3.2 and 4.0 pc. Upper panel: 
Evaporation + ejection. Lower panel: Evaporation. After birth 
the clusters spend a long time on the same spot along the dy- 
namical main sequence. For comparison the grey thick lines give 
the cooling track of the Standard Cluster (Fig. [§). 



100 
10 

1 

0.1 
0.01 
0.001 
0.0001 
100 
10 

1 

0.1 
0.01 
0.001 
0.0001 



Standard Cluster 
N=500 
N=2000 
N=3000 
N=4000 
N=5000 




Standard Cluster 
N=500 
N=2000 
N=3000 
N=4000 
N=5000 




0.01 



kT [M km 2 s" 2 ] 



Figure 13. Temperature-luminosity diagram for clusters of dif- 
ferent initial masses. Upper panel: Evaporation + ejection. Lower 
panel: Evaporation. Each cluster starts on the left and moves 
rapidly onto the common cooling sequence. For comparison the 
thick grey lines give the cooling track of the Standard Cluster 
(Fig- El- 



while the core is contracting and therefore the ansatz also 
does not hold. 

Furthermore, the findi ng that RhjR ude increases also 
contradicts the prediction o f lSpitzeij (| 19871 . p. 59), who found 
the ratio of Rh/Rtide to be a constant of about 3, i.e. the 
cluster to evolve self-similarly. Instead, this ratio is a func- 
tion of mass, and the value of 3 is reached just at the very 
end of a cluster's life-time. 

For the underlying set of investigated models, the con- 
stant value of the half-mass radius is a much better and 
more useful approximation. Therefore the half-mass radii of 
all clusters have been studied in detail and a constant has 
been fitted to R^ of each cluster between t cc and the disso- 
lution time, tdis- In Tab.[T]these fitted mean values, R^, for 
the half-mass radii of all model classes are shown. 



6.4 Tidal field 

The Galactocentric distance was set to 85 kpc for a number 
of models. Here the tidal force is much weaker and the tidal 
radius (eq. [8]) is about five times larger, which explains the 
following observed differences, since these are just due to a 
different scaling. 

Initially the cluster is almost a Standard Cluster with 
the same structure and relaxation time. It goes into core 
collapse on the same timescale, but the outer shells expand, 
due to the weaker tidal field, much further than those of 
a Standard Cluster (Fig. I17|l . This increases the dynamical 
time-scale of the cluster. After about 2 Gyr, the cluster has 
filled its tidal sphere and the half-mass radius then reaches 
a constant state. This shows that the tidal conditions have 
a strong influence on the equilibrium value of the half-mass 
radius. But taking a look at Fig. 1181 reveals that, after the 
equilibrium density distribution is established, the ratio of 
half-mass radius to tidal radius is the same for both clusters. 
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Figure 14. Temperature-luminosity diagram for the two clusters 
with masses of 16384 and 32768 Mq and the clusters with 5000 
Mq. Upper panel: Evaporation + ejection. Lower panel: Evapo- 
ration. Each cluster starts on the left and moves rapidly onto the 
common cooling sequence. For comparison the thick grey lines 
give the cooling track of the Standard Cluster (Fig. [6}. 
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Figure 15. The same as Fig. ll3l but with the cluster mass instead 
of the temperature on the x-axis. After core collapse all models 
are fully determined by their mass. Upper panel: Evaporation + 
ejection. Lower panel: Evaporation. For comparison the thick grey 
lines give the evolution of the Standard Cluster. 



Since encounters of stars are less frequent due to the 
smaller cluster density, and hence a lower temperature 
(Fig. I17p . the luminosity is much smaller than for a Standard 
Cluster (Fig. I19[) . The cooling tracks are therefore shifted 
within the T-L diagram. But in iV-body units the luminos- 
ity curves of the two cluster types would almost lie on top 
of each other, since after core collapse (or more precisely in 
this case: after both clusters have filled their tidal spheres) 
they just differ in terms of scaling. For example the total 
number of binary ejecters, which is independent of scaling, 
is therefore comparable to a Standard Cluster (also about 
10). 

The influence of a tidal field on cluster evolution can be 
observed in a direct comparison of a Standard Cluster with 
an isolated cluster of the same properties. Without a tidal 
field there is no tidal radius and therefore the criterion of 
escape has to be modified. Here, an escaper is defined as a 
star with positive energy, since this is the only way for a 



star to escape from an isolated cluster. The velocities of the 
escaping stars are not measured until they reach the 90% 
Lagrange radius. Otherwise the influence of the cluster po- 
tential on the stars' velocities would be too significant. Fur- 
thermore, scattering back to bound energies is very unlikely 
after reaching this radius. 

This difference in the escape criterion has some effect 
on the evolution in the T-L diagram (Fig. I19[) . Since the 
90% Lagrange radius is, during the first Gyr, much smaller 
than 2 Rude for the clusters at 85 kpc, the initial luminosity 
is above the level of those clusters because stars can fulfill 
the escape criterion more readily. But as cluster expansion 
goes on forever in an isolated cluster, due to the missing 
tidal field, escape is getting harder, since the 90% Lagrange 
radius also grows inexorably - and so does the half-mass 
radius. 

A tidal field therefore counteracts the cluster expan- 
sion, by removing the outer stars. An equilibrium value for 
the half-mass radius is established after core collapse, which 
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Figure 16. Ratio of half-mass radius to tidal radius for clusters 
with masses of 500 to 32768 Mq. The clusters asymptotically 
evolve onto a common sequence. Below 100 Mq fluctuations in 
Rh lead to a large scatter around this sequence. 



guarantees a balance between energy production in the core 
and energy loss through the tidal field. The absence of such 
a tidal field prevents the removal of weakly bound stars, so 
the cluster can expand further and further. 

Hence, after core collapse not only the number of stars 
left in the cluster are important but also the tidal conditions. 



6.5 Binaries 

Effects of primordial binaries on star cluste r evo- 
lution have been studied t horoughly l|Heggiel 1 19751 ; 
iMcMiUan. Hut fc Makinolll990l : lKroupalll995bT ). In this in- 
vestigation models with fbm = 0.95 are set up with 5% single 
stars and 95% centre of mass particles first. Thereafter the 
cm. particles are substituted with a random binary, follow- 
ing the binding energy distribution fiE (ea. I17[). 

Like in a Standard Cluster the initial temperature of sin- 
gle stars and cm. particles is about kT — 2M0km 2 s~ 2 . So, 
assuming that binaries with Ebin less than kT are disrupted 
immediately, setting up a binary fraction of fun = 0.95 with 
this binding energy distribution yields about 300 bound bi- 
nary systems initially in the cluster (Fig. I20[) . 

One major consequence of the presence of a large num- 
ber of primordial binaries is the well known effect of a 
damped core collapse. With a lot of energy sources available, 
the cluster has no need to go into deep collapse to produce 
energy fsee lOstrikeilll985l ). Furthermore, binary-binary in- 
teractions are much more efficient i n energy p roduction and 
therefore allow a lower core density l|Hutl 19961 ). Another way 
to put this is an analogon to a protostar which has a sig- 
nificant amount of primordial deuterium. Since deuterium 
burning is efficient at lower temperatures compared to hy- 
drogen burning, the whole star can achieve stability at a 
larger radius. 

This effect is accelerated by the fact that, since binary 
systems in single-mass clusters are subject to mass segre- 
gation, they sink to the cluster core within a few crossing 
times. So most binaries are burnt in the core within a few 
relaxation times at the very beginning of the computations 
and just a few remain (Fig. I20[) . 



100 



10 



2R tide 
LR90 




2000 4000 6000 8000 1 0000 1 2000 1 4000 1 6000 1 8000 
t [Myr] 

Figure 17. Above: Evolution of the Lagrangian radii for a cluster 
with Rq = 85 kpc. After core collapse the cluster is still expand- 
ing until it fills the tidal sphere at about 2000 Myr. The middle 
panel shows the mass evolution of the same cluster compared to 
a Standard Cluster, which is much slower due to the larger tidal 
radius and the increasing dynamical time-scale. In the bottom 
panel the temperature of the cluster is shown, which differs from 
the Standard Cluster after the latter has filled its tidal sphere, 
while the former continues expanding. 
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Figure 18. Ratio of half-mass radius to tidal radius for a clus- 
ter at Galactocentric distance of 85 kpc compared to a Standard 
Cluster at 8.5 kpc. 
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Figure 20. Upper panel: velocities of escaping stars as a function 
of time for a cluster with /j, in = 0.95. Below: number of binaries 
over the course of time. Most primordial binaries are burnt to 
prevent the core from collapsing and drive the halo expansion. 



the mass and the tidal conditions, but also by the binary 
fraction. 
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Figure 19. Mean temperature-luminosity diagram of an isolated 
cluster and one with a Galactocentric distance of 85 kpc. Up- 
per panel: Evaporation + ejection. Lower panel: Evaporation. For 
comparison the thick grey lines give the cooling track of the Stan- 
dard Cluster (Fig. [6}. 



The presence of primordial binaries also increases the 
relative number of ejecters, since there are many binaries 
to produce them already from the start (upper panel in 
Fig. 120 j l. Hence, the ejection luminosity is much larger than 
in a Standard Cluster while the evaporation is unaffected by 
this (Fig. I2ip . Nevertheless the limit between evaporaters 
and ejecters is the same as in a Standard Cluster. 

Furthermore, the first ejecter cannot be taken as a trig- 
ger for core collapse, contrary to the case of single-star clus- 
ters. But here also the standard criterion for core collapse of 
the first binary with a binding energy above a certain value 
naturally fails. 

Since most of the primordial binaries are disrupted in 
binary-binary encounters and just a few are ejected, much of 
their kinetic energy gets distributed among the other stars. 
This constant heat source leads to a more extended cluster 
and hence a different half- mass radius (see Tab. [I}. So after 
core collapse, a cluster may not be fully described by only 



7 CONCLUSIONS 

The investigations made here show that it is possible to 
define a dynamical temperature and a dynamical luminos- 
ity for star clusters. The dynamical temperature-luminosity 
diagram established with these two quantities gives a com- 
pletely new way of looking at cluster evolution and helps to 
gain valuable insights on the energy evolution of star clus- 
ters. 

The dynamical temperature is defined through the mean 
kinetic energy of the stars within a cluster, where binary sys- 
tems need to be replaced by their centre-of-mass particles. 
Through this correction the system can still be treated as 
an ideal gas, which enables to relate the temperature to the 
velocity dispersion of the stars within the cluster. In this 
way the dynamical temperature is directly correlated to the 
velocity distribution of the cluster stars, which is taken to 
be Maxwellian shaped. The latter implies that there is al- 
ways a high-velocity tail of stars with velocities above the 
escape velocity. When these unbound stars leave the cluster, 
the left-over stars will reestablish a velocity dispersion with 
a lower dynamical temperature within a relaxation time. A 
relation between dynamical temperature and the number of 
escaping stars is therefore expected, since there will always 
be a certain fraction of stars leaving the cluster. 

The dynamical luminosity is defined as the kinetic en- 
ergy of the stars leaving the cluster, in analogy to the energy 
of photons emitted by a star. Escaping stars show two dif- 
ferent origins, evaporation and ejection, which have to be 
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Figure 21. Mean temperature-luminosity diagram of the clus- 
ters with 95% primordial binaries. Upper panel: Evaporation + 
ejection. Lower panel: Evaporation. For comparison the grey lines 
give the evolution of the Standard Cluster (Fig. [6}. 



treated separately. The former is strongly correlated to the 
temperature, while the latter is due to binary interactions 
and hence a tracer for structural effects like core collapse 
or core oscillations. In fact, there is no ejection before core 
collapse and the occurrence of the first ejected star can be 
used to define the point in time when the core has reached its 
densest phase. The limit between evaporaters and ejecters 
was found to be universal among the investigated clusters 
at Eki n /(kT) = 5.1 ±0.2. These two families of escapers are 
furthermore the reason for splitting up the total luminosity 
into an evaporation luminosity and an ejection luminosity. 

In this way the energy evolution of a star cluster can 
be understood as follows: A cluster is generating energy in 
the core through binary burning which causes expansion of 
the whole cluster. Without a tidal field this leads to infinite 
expansion, since a single binary can generate large amounts 
of energy until it is ejected from the cluster, in which case 
the core has to contract to form a new one. Stars leaving 
the cluster have to have positive energy, which they can 
gain through two-body relaxation processes. This resulting 



evaporation of stars decreases the total energy of the clus- 
ter. Ejection mechanisms have to counteract this decrease 
of energy. Three- or more-body encounters in the core cause 
not only the ejection of stars with high velocities, but also 
the ejection of binaries. The latter therefore dissipate the 
accumulated binding energy and make the total energy of 
the cluster go to zero (see Fig. [5}. A tidal field is setting a 
limit to the expansion of the cluster, such that the cluster 
spills over the tidal radius. Evaporation is therefore ampli- 
fied, since escape is eased. 

All computed clusters show a tight relation between 
evaporation luminosity and dynamical temperature, while 
the ejection luminosity shows a much larger scatter due to 
the small number of ejecters, i.e. insufficient statistics. Be- 
fore core collapse all clusters move towards a dynamical main 
sequence, where they spend up to half of their total life-time. 
After core collapse they follow a common cooling sequence, 
similar to a cooling track of white dwarfs. The duration of 
the main sequence phase depends on the initial density dis- 
tribution of the cluster. For initially very extended clusters, 
i.e. with large initial Rh/Rude, core collapse takes a large 
fraction of the cluster's total life time, during which the 
temperature stays nearly constant, corresponding to a fixed 
position in the dynamical temperature-luminosity diagram. 

This leads to a large degeneracy in models for equal- 
mass clusters, since the initial conditions of a cluster on 
the cooling sequence cannot be traced back. For the evapo- 
ration luminosity the only deviations from this cooling se- 
quence are given by models with different tidal conditions, 
i.e. at different Galactic radii, while the ejection luminosity 
is additionally affected by the binary fraction. On the other 
hand this supports the hypothesis that after core collapse, 
the state of a single-mass system can be fully described by 
three quantities: the number of stars left in the cluster, the 
tidal conditions and the binary content. A theory of cluster 
evolution therefore has to focus on these parameters. 

Whether clusters with a stellar mass distribution would 
also form a dynamical main sequence still has to be inves- 
tigated. If this is the case, well observed clusters like the 
Pleiades may be placed within such a diagram and a unique 
evolutionary track may be assigned to them, which would 
give direct insights on the former and further evolution of 
those systems. 

Another important finding of this parameter-space 
study is the phenomenon of a constant half-mass radius for 
clusters in tidal fields. This holds for the time after core col- 
lapse, when a final density distribution is established within 
the tidal radius, in which the energy production in the core 
is balanced by the mass loss at the tidal radius. As shown in 
this work, the tidal field has a significant influence on this 
value, since clusters at larger Galactic radii show a much 
larger equilibrium value for the half-mass radius, while iso- 
lated clusters do not show this phenomenon at all. This can 
be expressed in terms of the ratio of half-mass radius to 
tidal radius, which evolves along a common sequence for all 
investigated clusters, depending only on mass. The ratio is 
increasing with time, as the half-mass radius is a constant 
and the tidal radius is decreasing due to ongoing mass loss. 
This implies that the energy production efficiency within the 
half-mass radius is increasing as the cluster loses mass. 

The Standard Cluster shows an equilibrium half-mass 
radius of about 2 pc in the given tidal field, and so do the 
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other models with the same initial mass, tidal conditions 
and binary fractions (Tab.[T]). An increasing initial mass in- 
creases the value of the equilibrium half-mass radius. The 
models with a higher mass than the Standard Cluster have 
a slightly declining half-mass radius at the end of their life- 
times, when their ratio of Rh/Rade has reached the common 
sequence. While the model with 16384 stars still shows a 
nearly constant Rh, the 32k model has a slightly increasing 
half-mass radius after core collapse until about half its disso- 
lution time, after which it starts decreasing slowly. In a first 
order approximation it still can be described by a constant 
but it is expected that this gets more and more inexact the 
larger the initial mass of the cluster is. 

This be haviour c ontra dicts the theory of self-similar 
evolution of ISpitzerl (|l987l . p. 59), who found the ratio 
Rh/Rude to be a constant of about 3, which seems to be 
a too rough estimate and just holds at the very end of a 
cluster's life-time. Furthermore, the simple cluster-evolution 
theory based on the ansatz that evaporation does not change 
the energy of a cluster and thus the half-mass radius scales 
with M 2 ()Binnev fe Tremainel Il987l . p. 525), can be com- 
pletely ruled out. This only holds for very extended clus- 
ters, for which the initial half-mass radius is larger than the 
equilibrium value, and for these clusters this ansatz holds 
only until core collapse because until then no ejection has 
occurred. In this given set of parameters the constant half- 
mass radius gives a much better approximation. Neverthe- 
less, this parameter-space study showed how similar single- 
mass clusters with a wide range of initial conditions evolve, 
once they have adjusted to the given tidal conditions, which 
is in most cases right after core collapse. 

The discovery of Russell and Hertzsprung was made 
possible through improvements in distance-determination 
techniques for stars, just like future developments like GAIA 
will increase the possibilities of measuring accurate peculiar 
velocities of stars in the Milky Way. Then, also a large sam- 
ple of open clusters will be surveyed in detail, which will 
give the opportunity to reconstruct internal quantities like 
the velocity dispersion very accurately and also to identify 
stars which are about to leave a cluster or have left it lately. 
This means that the actual measurement of the two quanti- 
ties, dynamical temperature and dynamical luminosity, will 
become feasible. 
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Table 1. Overview of all computed models, n gives the number of computed models of a particular kind, p(r) is the initial 

density profile (P: Plummer, K: King) and R^ is the fitted value for the half-mass radius after core collapse. Given errors 

are the standard deviations. The last column gives the fitted values for the limit between evaporaters and ejecters with 
fitting uncertainties. 

n N p(r) R h [pc] f bln R gal [kpc] t cc [Myr] t° el [Myr] t dis [Myr] r{ [pc] (^^.^ 

27 1000 P 0.8 8.5 165±20 10.1 2277±117 2.01±0.09 5.5lJ;° 

5 1000 K 0.8 8.5 162±37 10.1 2198±120 1.86±0.04 5. 4+^ 

5 1000 P 0.4 8.5 62±5 3.6 2039±85 1.96±0.12 5-7t\.l 

5 1000 P 1.6 8.5 432±50 28.7 2467±159 2.15±0.09 4.8+°* 

9 1000 P 2.4 8.5 758±47 52.6 2573±175 2.18±0.17 7.3+^ 
5 1000 P 3.2 8.5 1081±98 81.0 2592±173 2.09±0.09 5.9^2.2 

10 1000 P 4.0 8.5 1061±116 113.3 2489±94 2.03±0.13 4.3±og 

33 500 P 0.8 8.5 134±32 8.4 1662±174 1.88±0.12 6.7+.*;* 

8 2000 P 0.8 8.5 191±27 12.5 3195±150 2.13±0.05 ^t° % 

5 3000 P 0.8 8.5 241±11 14.2 3871±101 2.29±0.03 

4 4000 P 0.8 8.5 253±19 15.6 4563±316 2.38±0.14 4.3ljg 
3 5000 P 0.8 8.5 290±9 16.9 4943±180 2.48±0.13 4.4^^ 
1 16384 P 0.8 8.5 442 25.7 10077 2.99 3.8+^ 
1 32768 P 0.8 8.5 1093 33.3 15645 3.57 3.6^f 7 

5 1000 P 0.8 85 193±19 10.1 >20000 8.69±0.81 §.7t° 6 7 
5 1000 P 0.8 oo 146±31 10.1 >20000 (-^ 4.3+°5 



4 1000 P 0.8 0.95 8.5 7.8 2012±232 2.28±0.14 5.7 



Isolated clusters do not show a constant half-mass radius. 

The presence of primordial binaries prevents the core from collapsing. 



